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The color fields in a system composed by two static quarks and two static antiquarks are studied. 
In particular, we consider the four particles in the corners of a rectangle, and two possible alignment 
of the particles, one in which the quarks are at the same side of the rectangle, and the other where 
they are at opposite sides. We use a variational method, to probe not only the ground state but also 
the first excited state. This results permit us to observe and interpret the flux-tube recombination 
in the mesons to mesons and the tetraquark to mesons transitions, for both states. The results are 
compared with previous results for the static potential and the Casimir scaling predictions. 



I. INTRODUCTION 

Systems constituted by two quarks and two antiquarks 
are of extreme importance for strong interaction physics. 
Not only because they are a starting point for meson- 
meson scattering, but also because of the possible exis- 
tence of bound- states — tetraquarks, initially predicted 
by JafTe [1 . There are several observed resonances which 
are candidates to tetraquarks [2H1]. The most recent 
are the particles reported by the Belle collaboration 
[5j [6]. While it remains difficult to study tetraquarks 
in Lattice QCD with dynamical quarks [8] the static 
tetraquark potentials and flux tubes have been studied in 
precise lattice QCD simmulations [9HTT], The system of 
two quarks and two antiquarks in SU{2) was also studied 
using a variational method [12] . 

The Lattice QCD results for the static potential seem 
to confirm that the ground state potential is well de- 
scribed by a generalized flip-flop potential. This piece- 
wise potential was already proposed and used by different 
authors [13H2T]. 

as a device to eliminate the physically non existent 
long distance Van der Waals interactions arising from 
the naive potential models based on sum of two-body 
Casimir scaled potentials [22H27] . 

The lattice studies indicate that the ground state is 
well described by the generalized flip-flop potential 

V FF = mm(V I ,V II ,V T ). (1) 

This generalized potential differs from the firstly pro- 
posed flip-flop models because, besides the meson-meson 
domains already present in earlier models, it also includes 
a tetraquark domain. Vi and Vu, in Eq. 0, are the two 
possible two- meson potentials, given as sum of the two 
intrameson potentials. Explicitly, 

Vr = VM(|ri-r 3 |) + VM(|r2-r 4 |), (2) 

and 

V // = VM(|ri-r 4 |) + "M|r2-r 3 |). (3) 




Figure 1: The shape of the minimal string which links 
the four particles in the tetraquark. Note that when the 
diquark and the diantiquark are close the five segment 
structure collapses into a four segment one. 

where the intrameson potential Vm is well fitted in the 
static limit as a Cornell potential: V^(r) = C — ^ + err. 
The potential Vt, is the tetraquark potential and is given 
according to [10] and [9] by 

V T = C + aY J ^-^— +<?Lmin (4) 

where L min is the minimal distance linking the four par- 
ticles, as depicted in Fig. [T] 

However, the value of the static potential alone, is not 
sufficient to understand the confinement in tetraquarks. 
Confinement can be researched both with the measure- 
ment of the color fields and with the study of the colour 
wave functions of the tetraquark. Recently we computed 
the color-electric and color-magnetic fields generated by 
a static tetraquark [TT] . We confirmed that for quark 
and antiquark distances compatible with a tetraquark po- 
tential, tetraquark color flux tubes actually are created. 
Thus the mechanism of confinement in static tetraquarks 
is the localization of the color fields in fundamental flux 
tubes. 

But it remains important to clarify what are the color 
wave functions of static tetraquarks, and to observe the 
flip-flop of both the flux tubes and the wave functions 
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Figure 2: Wilson loop operator used by [9] and [10] to 
calculate the static potential of the tetraquark system. 



Q2 



(a) 



Q 4 Qa 

O A O 



ri 



(b) 



Q 2 



• o _ y • o 

Qi « r 2 Q3 Qi . r 2 Q3 



Figure 3: (a) Parallel alignment, (b) Antiparallel 
alignment. 



of static teraquarks. Here we apply the variational 
principle to observe whether the flip-flop transition oc- 
curs for the color flux tubes in a quar k- quark- ant iquar k- 
antiquark(QQQQ) system. 

In Section II, we discuss the two possible color singlets 
of the four-particle system, while in section III, we re- 
view the Wilson Loop of the tetraquark system. Then, 
in section IV, we make the synthesis of the two previous 
sections and develop a variational method to describe the 
QQQQ system, either in a tetraquark color state or in a 
two meson state. In section V, the method we use to com- 
pute the color fields in the lattice is described. In section 
VI, this method is used to compute the color fields in the 
system, for different arrangements of the four particles. 
In section VII our results are discussed and in section 
VIII the conclusions are presented. 



II. A TWO COMPONENT BASIS FOR THE 
TETRAQUARK COLOR WAVEFUNCTION 

Notice a tetraquark must be a color singlet due to con- 
finement and to gauge invariance, but there are two pos- 
sible color wave functions for a color singlet tetraquark. 
For instance, if the system is on the domain where 
Vff — Vj (domain I), we expect the color wavefunction 
of the ground state to be given by 



l^ZlQiQjQiQj) ■ 



(5) 



If the system is on domain //, we expect the ground state 
to be given by 



(6) 



Note that these two systems are not orthogonal as 

(i\n) = I 

When Vff = Vr, we expect that the two quarks form 
an antitriplet and the two ant iquar ks a triplet. So, the 



wavefunction is given by, 

\A) = Afeijk\QiQj)e k im\QiQr> 



(?) 



where the A stand for antisymmetric since this wave- 
function is antisymmetric for the exchange of two quarks 
or two antiquarks. Eq. ^ can be simplified by con- 
tracting the tensors, and by imposing the normalization 
(A\A) = 1, 
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We can also construct a color symmetric state as 



(\i)-\n))- 

(8) 



\S) 



2\/6 



(IQiQjQiQ^+lQiQjQjQi)) 



l(\i)+\n))- 

(9) 

Any color single state of two quarks and two antiquarks 
can be decomposed in two different sets of basis vectors, 
either \I) and \II) or \A) and \S). The second set of kets 
has the advantage of forming an orthonormal basis. \I) 
and \II) can be written as 



\I) = \ -o\S) 



1 
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\A), 



\H) 



\S) 



\A). 



(10) 



(11) 



Therefore, we can describe our state with a two com- 
ponent vector and the potential itself as a two by two 
matrix. The lowest eigenvalue of the matrix is corre- 
sponds to the static potential measured on the lattice, 
with the corresponding eigenvector being \II) or \A) 
depending on the domain considered. Since the potential 
is hermitian, the eigenvector of the excited state of the 
potential is orthogonal to the one of the ground state. 
Therefore, these eigenvectors must be 



1 



\S) 



\A), 



(12) 



3 




Figure 4: Elements of the Wilson loop matrix in the 
parallel alignment case. 



Figure 5: Elements of the Wilson loop matrix in the 
antiparallel alignment case. 



for the domain I, and 



(13) 



for the domain II. In the tetraquark domain the excited 
state eigenvector must be \S), defined in Eq. 



III. THE TETRAQUARK WILSON LOOP 

The static potential for the tetraquark has been stud- 
ied in the lattice by [10J and |9J. The Wilson loop oper- 
ator for the tetraquark system, illustrated in Fig. [2j is 
given by 



W 4Q = -Tt[M 1 R 12 M 2 L 12 ] . 



with 



L 



tdii _ * c c t>33 r>kk' 



(14) 



(15) 
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This Wilson Loop has the quantum numbers of color 
singlet system where the two quarks form an antitriplet 
and the two antiquarks form a triplet. 

This lattice studies, indicate that the static potential of 
the two-quark and two-antiquark system is a generalized 
flip-flop potential, 



Vff = min(Vr, Vm 1 m 2 , Vm 3 m 4 ) : 



(16) 



where Vm 1 m 2 an d Vm 3 m 4 are the two possible two-meson 
potentials, given by the sum of two independent intra- 
meson potentials Vm 1 m 2 — Vm x + Vm 2 , an d is the 
tetraquark potential which corresponds to the sector 
where the four particles are confined, linked by a single 
fundamental string. 

For the special case where the four particles form as 
rectangle as in Fig. [TJ 

Lmin is given by — y/3vi -\-V2 
for r 2 > ^ [H [29]. If we neglect the Coulomb part of 
Vff, the flip-flop potential with the linear potentials only 
produces a tetraquark domain for distances r 2 ^ y/Sri • 



IV. VARIATIONAL METHOD FOR THE QQQQ 
SYSTEM 

Now we extend the Wilson loop in order to obtain 
the true ground state of the tetraquark system and also 
to obtain the first excited state. To achieve this, we 
note that the Wilson Loop operator can be written as 
a correlation of a certain operator at different times 
W(t) = (O(t)dt(o)). This can be generalized by con- 
sidering instead a base of operators Oi. So this way, our 
Wilson loop becomes a matrix Wij = (OiOj). This can 
be used not only to improve the ground state overlap but 
also to obtain more energy levels of the system. The en- 
ergy levels are the solution of the generalized eigensystem 



(W ij (t))ci(t)=w n (W ij (0))ci(t). 



(17) 



We consider again the case where the four particles 
form a rectangles and also two different alignments of the 
QQQQ system. A parallel one, where the two quarks are 
on the same side of the rectangle and an antiparallel one, 
where the quarks are on opposite corners of the rectangle, 
see Fig. [3] 

For both cases we use a basis of two operators, in- 
spired in Ref. [30J. In the parallel case the opera- 
tors are the tetraquark operator and a two me- 
son operator. This gives a Wilson loop matrix where 
the diagonal elements are: the tetraquark Wilson Loop 
W 4 q = (04q(*)04q(O)) (Fig. [i]), the correlation of two 
Wilson loops (one per meson composed of a quark and 
an antiquark) and the off-diagonal elements correspond- 
ing to the transition between the two states. The four 
matrix elements, each one corresponding to a different 
Wilson loop, are depicted in Fig. [4j 

For the antiparallel alignment, the operators we uti- 
lize are the two meson- meson operators, giving the four 
matrix elements given on Fig. |5j 



V. CHROMO-FIELDS COMPUTATION 

We compute the color electric and the color magnetic 
fields, utilizing the correlators of the plaquettes P^ v and 
the Wilson loops W n . We define the plaquettes as P^ v = 
1 - \Tr\U^*)U v (s + /x)C/t( s + „)tft( B )]. 
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With this definition, the chromo- fields are given by 



(Ef) = (P 0i ) 



(W n P 0i ) 
(W) 



/b2) _ (WP jk ) _ 

K (w n ) ™' 



(18) 
(19) 



with the indices j and k complementing index i. The La- 
grangian and Energy densities are given by C = \(E 2 — 
B 2 ) and U = \{E 2 + B 2 ). The Wilson loop operator for 
each state is given by W n (t) = c % n (t)WijC J n (t), in confor- 
mity with Eq. ( fl7| ). 

The plaquettes are placed at t = T/2, where T is the 
temporal extension of the Wilson Loops. We calculate 
the correlators for different values of T and fit the results 
from T = T min to T = T max to a constant, and calculate 
the x 2 /d.o.f.. We find that for T m ^ n = 4 and T max = 16, 
the values of the x 2 /d.o.f. are generaly acceptable, while 
still giving a clear signal. For this we consider the value 
of the fields as the result of the fit to a constant of the 
correlators from T = 4toT = 16. 



VI. RESULTS FOR THE COLOR FIELDS 

In this work, we only consider the geometry of the 
four particles in a plane forming a rectangle with two 
distinct kinds of alignment, see Fig. [3J When the two 
quarks/antiquarks are on the same side of the rectangle, 
this corresponds to a parallel alignment; if they are in 
opposite sides, we have an antiparallel alignment. The 
parallel alignment is useful in the visualization of the 
transition between the tetraquark and the two meson 
states, while the antiparallel one is be used to observe 
the transition between the two meson- meson states. 

The results presented in this work are obtained us- 
ing 1121 quenched lattice QCD configurations with di- 
mension 24 3 x 48 and (3 = 6.2 with lattice spacing, a, 
a = 0.07261(85) fm or a' 1 = 2718(32) MeV The configu- 
rations were generated with GPUs using a combination of 
Cabbibo-Marinari, pseudo-heatbath and over-relaxation 
algorithms, [3T| [32], APE smearing [33] and Hyper cubic 
blocking [34] are used to improve the signal to noise ratio. 

In Figs. [6] and [7| we show the results of the La- 
grangian density for the antiparallel alignment for the 
ground state, n = 0, and the first excited state, n = 1, 
respectively, with different values of T\ and 7*2. Note 
that, in this geometry the system is symmetric for the 
exchange of T\ and 7*2. When r\ = 7*2, the system forms 
a square-symmetric structure for both the ground and 
first excited states, thus the two states have to corre- 
spond to the symmetric and antisymmetric color wave- 
functions. Inspecting the composition of the variational 



Wilson Loop for the two states, we are able to conclude 
that the ground state corresponds in this geometry to 
the color symmetric wavefunction, with the color anti- 
symmetric corresponding to the first excited state. This 
is the reverse of what happens in the tetraquark sector, 
where the ground state is color antisymmetric. When 
7*1 < r2, the groundstate results show a pure two-meson 
ground state. In what concerns the first excited state, it 
must have a color wavefunction orthogonal to the one of 
the ground state. We observe a first excited state differ- 
ent from the other two- meson state, where the confining 
string seems to be linking all the four particles together. 

In Figs. |8]and|9j we depict the plot of the Lagrangian 
density of the parallel alignment, in the ground and first 
excited states respectively. These results are comple- 
mented with the Figs. 
x = and y 



m and 11 



0, for 7*1 = 6 and 7*2 



with cuts cuts for 
8. 



We now analyse the ground states profiles. In upper 
graphics of both Fig. [lO]and[TT] we see the results for the 
cuts of the Lagrangian density (with 7*1 = 6 and 7*1 = 8 
respectively), for x = (on the left side) and y = (on 
the right side). The first case corresponds to the center 
of the predictable diquark-diantiquark flux tube formed 
in the tetraquark domain, while the second corresponds 
to a transversal cut of the same flux tube and also of the 
two mesonic flux tubes, whose the centers should be in 
— 7*i/2 and 7*2/2. We observe the transition between the 
tetraquark and the two meson states. 

More interesting and not so easily interpretable is the 
first excited state. We can with a fair amount of certainty 
claim that this first excited state is due to the recombi- 
nation of the flux tubes and not due to the flux tube 
excitations as in j35j|36]. For the distances studied, we 
always seem to have a fully confined system, with all the 



four particles connected by the flux tube (Figs. [Tl [9^, 10 



and 11) contrarily to what we can naively be expect for 
the meson to meson transition. The first excited state 
clearly is not the other two meson state, but a different 
color state. Neither, is it a vibrational excitation of a 
flux tube. 

Indeed the first excited state should be orthogonal in 
the color space to the ground state. Since both the 
ground state and the first excited state are eigenvectors of 
the static potential matrix since the phenomenon which 
is happening is the recombination of the flux tube, and 
since this matrix is hermitian, the two vectors have to 
be orthogonal to each other. Consequently, we expect 
that 1 7) to be the color vector of the first excited state, 
when the ground state has the color of the meson system 
I J). Similarly, we expect that the first excited state to be 
given by the symmetric color state \S), when the ground 
state is the antisymmetric one \A). This happens when 
we are in the tetraquark domain. 
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L for r 1 = 8 and r 2 = 8 




2.54e-03 
i:80e-O3 


8 


S.OOe-04 


6 


8.00e-04 


4 


7.00e-04 


2 


6.00*44 




5.00e-04 




4.00e-04 


-2 


3.00e-04 


-4 


2.00*44 




1.00e-04 


-6 


040*440 


-8 


2.53*43 
l2>Oe-03 


8 


9.00*44 


6 


8.00e-04 


4 


740*44 


2 


6.00*44 




5.00e-04 




440*44 


-2 


3.00*44 


-4 


2.00e-04 




1.00e-04 


-6 


0.00*400 


-8 






2.53e-03 


8 




= 


2.53e-03 




9.00e-04 


6 






9.00e-04 




8.00e-04 


4 






S.00e-04 




7j00e-04 


2 






7.00*44 




6.00*44 








6.00*44 




Sj00e-04 


>.o 






5.00e-04 




4.00e-04 


-2 






4.00e-04 




3.00e-04 


-4 






3.00e-04 




2.00e-04 






I 


2.00*44 




1.00e-04 


-6 






1.00e-04 


: 


0.00e+00 


-8 






0.00e+00 



L for r, =8 and r 2 = 10 



-8 -6 -4 -2 2 4 



L for r 1 = 10 and r 2 = 10 






8 




9.00*44 


6 




8.00e-04 


4 




740*44 


2 




6.00e-04 






5.00e-04 


>.o 




4j00e-04 


-2 




3.00e-04 


-4 




2.00e-04 






1.00e-04 


-6 


: 


0.00*400 


-8 




Figure 6: Lagrangian density for the ground state of the antiparallel alignment. 



L for r- = 6 and r. 




L for r 1 = 8 and r 2 = 8 
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L for r, =8 and r 2 = 10 
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Figure 7: Lagrangian density for the first excited state of the antiparallel alignment. 



VII. DISCUSSION 



As we discussed before, the results seem to agree 
with the previous ones, obtained for the static potential 
[9j [10] , which support the generalized flip-flop picture for 
the ground state of two quarks and two antiquarks sys- 
tem. Namely, we observe the formation of the tetraquark 



string and of the two mesonic strings in the domains we 
would expect that to happen. 

For the first excited states, the situation is not as clear, 
and so to better understand these states we calculate the 
Casimir scaling factors for the different color wavefunc- 
tions in Table |TJ This shall provide us with qualitative 
insight to what kind of interaction we expect in the sys- 
tem. 
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Figure 8: Lagrangian density for the ground state of the parallel alignment. 
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Figure 9: Lagrangian density for the first excited state of the parallel alignment. 



In Table [T] we see that Casimir Scalling predicts a 
repulsive quark- quark and antiquark-antiquark interac- 
tion and an attractive quark-antiquark interaction for the 
state | S). As can be seen in Figs. [9, 10 and 11 , the results 
qualitatively agree with this prediction, due to a suppres- 
sion of the flux-tube between the two quarks and the two 
antiquarks. For the state |J), the prediction is the repul- 
sion between the particles that form the two mesons in 
the ground states, while all other interactions are attrac- 
tive. Again, the results do not qualitatively contradict 
this possibility, as shown in Fig. [9j for (ri,V2) = (6,8) 



and (r*i, T2) = (8, 10), where we again see a flux tube sup- 
pression between the particles that form mesons in the 
ground state. 



So, the results for the first excited state seem to agree 
with the hypothesis that the first excited state is orthog- 
onal to the ground state in the color space, with \I) and 
\II) and \S) being the excited states, where |J), \II) and 
I A) respectively are the ground states. 
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Figure 10: Cuts of the Lagrangian density for x = and for y = 0, with n = 6, both for the ground state and for 

the first excited state. 



VIII. CONCLUSION 

In this work, we use a variational method to com- 
pute the chromo-fields of the system composed of two 
quarks and two antiquarks. With this method we can 
not only observe the region where the tetraquark state 
is the ground one, but also the region where the ground 
state is composed of two mesons, as well as the transi- 
tions between this two regimes and the onw between the 
two possible meson- meson ground states. We are also 
able to observe, for the first time, the first excited state 
of this system. 

For the ground state the results improve our previous 
work [11] where only the tetraquark operator was used. 
There we had difficulties with the measurement of the 
ground state outside of the tetraquark region due to the 
low overlap of the used operator with the true ground 
state. Here, we are able to overcome that difficulty by 
using a variational basis. The results are similar, in the 
tetraquark region, to those obtained there, while giving 
the expected transition to the meson-meson behavior out- 



side that region. 

The results for the first excited state are not, at the 
first view, as understandable as the ones for the ground 
state. We note however that since the ground state is well 
bellow the first gluonic excitation of the string, the only 
way by which we can explain this excitation is the flux 
tube recombination. This way we know that the first ex- 
cited state has a color eigenfunction which is orthogonal 
to the ground state, be it a tetraquark or a two-meson 
state. By using, the Casimir factors, we compare the 
predictions of the aforementioned hypothesis with our 
results of the Lagrangian density for this state. We con- 
clude that both results are compatible. 

Note that this first excited state should be as impor- 
tant as the ground state of this system. Think for in- 
stance in the decay of a tetraquark into a two meson 
system. Since the color structure of the initial and final 
states are not the same, we can not consider the potential 
of the system to be a scalar in the color space. Instead 
it should be given by a two by two matrix as explainded 
above. To reconstruct this matrix we need to know not 
only the potential of the ground state but also it's color 
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Figure 11: Cuts of the Lagrangian density for x — and for y = 0, for n = 8, both for the ground state and for the 

first excited state. 
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composition as well as the potential of the first excited 
state. 



This study should be complemented in the future by a 
detailed study of the potential and color composition of 
the first excited state in the lattice. 
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